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Description 

INCORPORATION OF UNCERTAINTY 
INFORMATION IN MODELING A 
CHARACTERISTIC OF A DEVICE 

BACKGROUND OF INVENTION 

[0001] 1. Technical Field 

[0002] The present invention relates to a method and model for incorporating 
uncertainty information when modeling a characteristic of a design. 

[0003] 2. Related Art 

[0004] An integrated circuit is subject to process variations that occur during 

manufacture of integrated circuits. Circuit designers currently assess the 
effect these process variations on performance characteristics of 
integrated circuits by methods which do not yield sufficiently accurate 
assessments. Thus, there is a need for a method that assesses the effect 
of process variations on performance characteristics of integrated circuits 
more accurately than is accomplished in the related art. 

SUMMARY OF INVENTION 

[0005] The present invention provides a method for modeling a characteristic C 
that is distributed within a domain, said method comprising: 

[0006] 

providing a base equation expressing the characteristic C as a function f 



of a variable V through use of N+1 parameters C 0 , C v .... C N , said base 

equation being of the form C = f(C 0 , C 1 C N , V), said N being at least 1, 

said parameters C 0 , C r .... C N being subject to uncertainty; 

[0007] providing a probability density function (PDF) for describing the probability 
of occurrence of C 0 in accordance with said uncertainty; and 

[0008] providing subsidiary equations expressing C 1 C N in terms of C 0 . 

[0009] The present invention provides a method for modeling a characteristic C 
that is distributed within a domain, said characteristic C having J 
subcharacteristics S.,, S 2 S Jf said method comprising: 

[0010] providing a combination equation that expresses C as a function F of the J 
subcharacteristics, said J being at least 2; 

[0011] providing base equations expressing Sj as a function fj of a variable V 

through use of N+1 parameters S j0 , S^, S jN , said base equations being 

of the form Sj = fj(S j0 , S j1 S jN , V), said N being at least 1, said 

parameters S J0 , S j1f S jN being subject to uncertainty, said j having 
values of 1, 2, .... J; 

[0012] providing at least one probability density function (PDF) from PDF V PDF 2 , 
.... PDF J( said PDF n describing the probability of occurrence of S n0 in 

accordance with said uncertainty for n=1 , 2, J, said at least one PDF 
including PDF.,; 

[0013] for each PDF n not provided: providing an auxiliary equation E n expressing 
S nn in terms of S 1fl ; and 



[0014] providing subsidiary equations expressing S j1 S jN in terms of S j0 for 

each subcharacteristic Sj of the J subcharacteristics. 

[0015] The present invention provides a computer program product, comprising a 
computer usable medium having a computer readable program code 
embodied therein, said computer readable program code adapted to be 
executed on a processor for implementing a method for modeling a 
characteristic C that is distributed within a domain, said method 
comprising: 

[0016] providing a base equation expressing the characteristic C as a function f 
of a variable V through use of N+1 parameters C 0 , C v C N , said base 
equation being of the form C = f(C 0 , C v C N , V), said N being at least 1, 
said parameters C 0 , C v C N being subject to uncertainty; 

[0017] providing a probability density function (PDF) for describing the probability 
of occurrence of C 0 in accordance with said uncertainty; and 

[0018] providing subsidiary equations expressing C v C N in terms of C 0 . 

[0019] The present invention provides a computer program product, comprising a 
computer usable medium having a computer readable program code 
embodied therein, said computer readable program code adapted to be 
executed on a processor for implementing a method for modeling a 
characteristic C that is distributed within a domain, said characteristic C 
having J subcharacteristics S 1t S 2 Sj, said method comprising: 



[0020] 

providing a combination equation that expresses C as a function F of the J 



subcharacteristics, said J being at least 2; 

[0021] providing base equations expressing Sj as a function f. of a variable V 

through use of N+1 parameters S j0 , S j1 S jN , said base equations being 

of the form S. = fj(S j0> S jN , V), said N being at least 1, said 

parameters S j0 , S j1 S jN being subject to uncertainty, said j having 

values of 1, 2, .... J; 

[0022] providing at least one probability density function (PDF) from PDF V PDF 2 , 
.... PDFj, said PDF n describing the probability of occurrence of S n0 in 

accordance with said uncertainty for n=1 , 2, .... J, said at least one PDF 
including PDF.,; 

[0023] for each PDF n not provided: providing an auxiliary equation E n expressing 
S n0 in terms of S 10 ; and 

[0024] providing subsidiary equations expressing S^, S jN in terms of S j0 for 
each subcharacteristic Sj of the J subcharacteristics. 

[0025] The present invention provides a model, comprising: 

[0026] a base equation expressing a characteristic C as a function f of a variable 
V through use of N+1 parameters C 0> C r .... C N , said base equation being 

of the form C = f(C 0 , C 1 C N , V), said N being at least 1 , said 

parameters C 0> C v .... C N being subject to uncertainty, said characteristic 
C being distributed within a domain; 

[0027] a probability density function (PDF) for describing the probability of 
occurrence of C 0 in accordance with said uncertainty; and 



[0028] subsidiary equations expressing C v .... C N in terms of C Q . 
[0029] The present invention provides a model, comprising: 

[0030] a combination equation that expresses C as a function F of the J 

subcharacteristics, said J being at least 2, said characteristic C being 
distributed within a domain, said characteristic C having J 
subcharacteristics S 1( S 2 , Sj; 

[0031] base equations expressing Sj as a function fj of a variable V through use 

of N+1 parameters S j0 , S j1 S jN , said base equations being of the form 

Sj = f^SjQ, S^, S jNI , V), said N being at least 1, said parameters S j0 , S^, 
.... S jN being subject to uncertainty, said j having values of 1 , 2, J; 

[0032] at least one probability density function (PDF) from PDF.,, PDF 2> .... PDFj, 
said PDF n describing the probability of occurrence of S n0 in accordance 

with said uncertainty for n=1 , 2 J, said at least one PDF including 

PDF.,; 

[0033] for each PDF n that does not exist: an auxiliary equation E n expressing S n0 
in terms of S 10 ; and 

[0034] subsidiary equations expressing S j1 S jN in terms of S j0 for each 

subcharacteristic Sj of the J subcharacteristics. 

[0035] The present invention provides a method that assesses the effect of 
process variations on performance characteristics of integrated circuits 
more accurately than is accomplished in the related art. 

BRIEF DESCRIPTION OF DRAWINGS 



[0036] FIG. 1 depicts a nominal curve of capacitance versus applied voltage of a 
varactor fitted to a fifth order polynomial in the applied voltage, in 
accordance with embodiments of the present invention. 

[0037] FIG. 2 depicts the effect of process variations on the nominal curve of FIG. 
1 based on six independent degrees of statistical freedom, in accordance 
with embodiments of the present invention. 

[0038] FIGS. 3A-3E depict curves relating parameters used for modeling the 
capacitance of a varactor with each other, in accordance with 
embodiments of the present invention. 

[0039] FIG. 4 depicts the effect of process variations on the nominal curve of FIG. 
1 based on one independent degree of statistical freedom, in accordance 
with embodiments of the present invention. 

[0040] FIG. 5 is a flow chart for providing a model for a characteristic C 

distributed in a domain wherein the characteristic C is not expressed in 
terms of subcharacteristics of C, in accordance with embodiments of the 
present invention. 

[0041] FIG. 6 is a flow chart for randomly sampling C based on the model 
provided in FIG. 5, in accordance with embodiments of the present 
invention. 

[0042] FIG. 7 is a flow chart for determining a performance parameter of a design 
based on the model of C in FIG. 5 and on the method for sampling C in 
FIG. 6, in accordance with embodiments of the present invention. 



[0043] FIG. 8 is a flow chart for providing a model for a characteristic C 

distributed in a domain wherein the characteristic C is expressed in terms 
of subcharacteristics of C, in accordance with embodiments of the present 
invention. 

[0044] FIG. 9 is a flow chart for randomly sampling C based on the model 
provided in FIG. 8, in accordance with embodiments of the present 
invention. 

[0045] FIG. 10 is a flow chart for determining a performance parameter of a 

design based on the model of C in FIG. 8 and on the method for sampling 
C in FIG. 9, in accordance with embodiments of the present invention. 

[0046] FIG. 1 1 illustrates a computer system used for implementing one or more 
algorithms that model a characteristic of a design, in accordance with 
embodiments of the present invention. 

DETAILED DESCRIPTION 

[0047] 

The present invention relates to the design of a physical system 
comprising devices therein such that characteristics of the devices are 
subject to uncertainty. For example, a physical system may comprise 
integrated circuits on a semiconductor wafer, wherein characteristics (e.g., 
capacitance) of devices (e.g., virtual capacitors) on the wafer may be 
subject to uncertainty due to process variations that occur during 
manufacturing of the integrated circuits on the semiconductor wafer. To 
confirm that such a design is acceptable, the designer (e.g., the integrated 
circuit designer) may assess the design by computerized simulation of the 



system for determining whether performance characteristics of the system 
are within acceptable limits. For integrated circuit design, such 
performance characteristics may include, inter alia, signal transmission 
speed, signal switching speed, operating temperature, power output, 
power conversion efficiency, the logical correctness of a logic circuit, etc. 
The computerized simulation takes into account the uncertainty, such as 
uncertainty due to process variations, in the device characteristics by the 
methodology of the present invention. The uncertainties are modeled by 
one or more probability density functions, such as a normal probability 
distribution characterized by a mean and a standard deviation. In addition 
to capacitance, other examples of device characteristics in electrical 
systems include, inter alia, the resistance of a metal line, the inductance of 
an inductor, and the threshold voltage of a transistor. 

[0048] 

The subsequent description of the present invention is divided into four 
segments. The first segment illustrates the present invention by discussing 
the modeling of the characteristic C of the capacitance of a varactor, 
which is a capacitor whose capacitance varies with the applied voltage. 
The second segment formulates the present invention generally by 
discussing the modeling of a characteristic C that is distributed within a 
domain, wherein C is not expressed in terms of any subcharacteristics of 
C, as embodied in Equation (2) discussed infra. The third segment 
formulates the present invention generally by discussing the modeling of a 
characteristic C that is distributed within a domain, wherein C is expressed 
in terms of a plurality of subcharacteristics of C, as embodied in Equation 



(3) discussed infra. The fourth segment describes a computer system that 
may be utilized to implement algorithms of the present invention. 

[0049] Modeling the Capacitance of a Varactor 

[0050] A particular example of a device in an integrated circuit is a varactor built 
in a semiconductor via by implanting and diffusing dopant impurities to 
create a p-n junction. Such a varactor may be described by Equation (1): 

[0051] C = C</(1-V/V B ) m (Equation"!) 

[0052] where C is the capacitance of the varactor, V B is the 'built-in' voltage 
(V B >0), V is the applied voltage at the varactor (V< 0), and m is the 
"grading coefficient" (m>0). Note that C = C 0 if V = V 0 =0. The use of the 
variable "C" for capacitance in Equation (1) and similar equations should 
not be confused with the use of "C" as representing a characteristic 
generally in the description of the present invention herein; e.g., in 
Equation (6) described infra. The context of the pertinent description in 
which "C" appears makes it clear when "C" represents capacitance and 
when "C" represents a characteristic generally. 

[0053] Tne varac tor described by Equation (1) is a "virtual capacitor" as 
contrasted with a physical capacitor that has two conducting plates 
separated by a dielectric medium such as air. A virtual capacitor has 
virtual capacitor plates with an effective area and an effective perimeter. 
For the case of virtual capacitance caused by dopant concentrations in 
semiconductors, the effective plate area represents the region of the p-n 
junction having uniform depth, and the effective plate perimeter represents 



the regions of the diffused p-n junction which is not of uniform depth, 
typically on the edges of the diffused region. A real or virtual capacitor 
may be characterized by an area capacitance and a perimeter 
capacitance. Area capacitance is defined as capacitance per unit area of a 
physical capacitor plate if the capacitor is a physical capacitor, or the 
capacitance per unit effective area of a virtual capacitor plate if the 
capacitor is a virtual capacitor. Perimeter capacitance is defined as 
capacitance per unit perimeter of a physical capacitor plate if the capacitor 
is a physical capacitor, or the capacitance per unit effective perimeter of a 
virtual capacitor plate if the capacitor is a virtual capacitor. Generally, the 
scope of the present invention includes both physical devices and virtual 
devices. 

[0054] In Equation (1), the parameters C 0 , V B , and m are subject to uncertainty, 
such as due to process variations, and are modeled statistically in 
accordance with the present invention. Note that Equation (1) has a 
physical basis for describing the capacitance C of a varactor and may 
accurately describe a varactor in many situations. In some situation, 
however, an equation with more adjustable parameters than the three 
parameters (C 0 , V 0 , and m) in Equation (1) may be useful for fitting 

measured capacitances. Accordingly, Equation (2) describes the 
capacitance C of a particularly useful varactor, namely a hyperabrupt 
varactor, using a fifth order polynomial in V with fitting six coefficients of 
C 0 , C v C 2 , C 3 , C 4 , and C 5 . 

[0055] C = C 0 (1 + C,V + C 2 V 2 + C 3 V 3 + C 4 V 4 + C 5 V 5 ) (Equation 2) 



[0056] Note that C = C 0 if V = V p =0. As in Equation (1), V is the applied voltage. 
Unlike Equation (1), the fifth order polynomial in V of Equation (2) has no 
physical basis for describing the hyperabrupt varactor and is utilized 
because the fifth order polynomial in Equation (2) can be accurately fitted 
to measured total capacitances of hyperabrupt varactors as a function of 
applied voltage, as illustrated in FIG. 1 in accordance with embodiments of 
the present invention. In the testing associated with the measurements of 
FIG. 1, the total capacitances of 5 devices on each of 12 dies were 
measured from which the data shown in FIG. 1 are drawn. This measured 
data was used to extract area and perimeter capacitances for each 
device. Each die of the 12 dies had five devices and one device was 
selected from each die for the data shown in FIG. 1 . All 5 devices were 
located close to one another (i.e., within about 3000 microns) on the same 
die. The anode to cathode voltage of the varactor is swept from slightly 
forward bias to relatively large reverse bias but less than the onset of 
avalanche multiplication. 

[0057] FIG. 1 shows that the nominal curve 10 fitted in accordance with Equation 
(2) accurately models the measured capacitances. Although the nominal 
curve 1 0 of FIG. 1 , which is fitted in accordance with Equation (2), 
accurately models the measured capacitances, the nominal curve 10 does 
not account for uncertainties in the capacitance C due to process 
variations. 

[0058] 

To account for the process variations, the consequent uncertainties in C 
may be simulated by attributing uncertainties to the parameters C 0 , C v C 2 , 



C 3 , C 4 , and C 5 of Equation (2). Accordingly, FIG. 2 illustrates capacitance 
versus voltage curves obtained by selecting C 0 , C v C 2 , C 3 , C 4 , and C 5 as 
being randomly distributed about a nominal fit to Equation (2) (e.g., the 
nominal fit 10 of FIG. 1 described supra), in accordance with embodiments 
of the present invention. The parameters C 0 , C v C 2 , C 3 , C 4 , and C 5 were 
randomly sampled independently and the resultant values of C 0 , C r C 2 , 
C 3 , C 4 , and C 5 were substituted into Equation 2 to generate the curves of 
FIG. 2. The curve 20 represents a nominal fit to Equation (2), whereas 
curves 21-29 were generated using sampling constraints. In particular, 
each coefficient of the set of coefficients (i.e., C 0 , C v C 2 , C 3 , C 4 , and C 5 ) 

was independently randomly selected from a normal distribution 
characterized by a mean equal to the nominal value of the coefficient and 
a standard deviation a such that 3a is 10% of the nominal value of the 
coefficient (e.g., if the nominal value is 3 then a is 0.1). Thus, each curve 
of curves 21-29 has a unique set of coefficients randomly selected as 
described supra, such that the voltage was varied for the unique set of 
coefficients. 

[0059] 

Although the nominal curve 20 in FIG. 2 is physically realistic, curves 21- 
29 do not accurately model the uncertainties in the capacitance C and 
include highly non-physical capacitances including negative capacitance 
values. Thus, incorporating random process variations into the model of 
Equation (2) by allowing random, independent variations of the 
parameters C 0 , C.,, C 2 , C 3 , C 4 , and C 5 does not result in a physically 

realistic description of the effect of random process variations on such a 



varactor. This non-physical result might be expected, since Equation (2) 
has no particular physical derivation; i.e., the polynomial form of Equation 
(2) was chosen for its utility in fitting the experimental curves. 
Nonetheless, the present invention modifies the preceding methodology in 
a way that generates physically realistic description of the effect of random 
process variations while retaining the form of Equation (2), as will be 
described next. 

[0060] The inventors of the present invention have determined that the non- 
physical nature of the curves 21-29 in FIG. 2 is caused by randomly 
varying the parameters C 0 , C r C 2 , C 3 , C 4 , and C 5 independently, which is 

equivalent to having six degrees of statistical freedom. The six 
independent degrees of statistical freedom is incorrect because to have a 
physically realistic result there should be only one degree of statistical 
freedom. The basis of having only one degree of statistical freedom is 
supported by FIGS. 3A-3E which shows five correlating relationships 
among the six parameters C 0 , C 1f C 2 , C 3 , C 4 , and C 5 , in accordance with 
embodiments of the present invention. In other words, if C 0 is permitted to 

vary randomly in conjunction with Equation (2), then the other parameters 
C r C 2 , C 3 , C 4 , and C 5 are determined by being correlated with C 0 , as 

shown for the area capacitance coefficients plotted in FIGS. 3A-3E. 

[0061] 

In the tests used to infer the plotted data curves in FIGS. 3A-3E, the total 
capacitances of 5 devices on each of 42 dies were measured from which 
the area capacitance coefficients shown in FIGS. 3A-3E are drawn. This 
measured data was used to extract the area and perimeter capacitances. 



All 5 devices were located close to one another (i.e., within about 3000 
microns) on the same die. All 42 dies were used to generate the area 
capacitance coefficient data shown in FIG. 1. Forty-two sets of polynomial 
coefficients were determined from these 42 die to model the area 
capacitance and another 42 sets of polynomial coefficients were 
determined from these 42 die to model the perimeter capacitance. 
Regression analyses of the polynomial coefficients used to model the area 
capacitance are shown in FIGS. 3A-3E. 

[0062] FIGS. 3A-3E each show a data curve and a curve representing a 

quadratic fit to the data. FIG. 3A depicts C, versus C 0 , FIG. 3B depicts C 2 
versus C v FIG. 3C depicts C 3 versus C 2 , FIG. 3D depicts C 4 versus C 3> 
and FIG. 3E depicts C 5 versus C 4 . Thus, the correlations shown in 
FIGS.3A-3E are equivalent to C.,, C 2 , C 3 , C 4 , and C 5 being each correlated 
with C 0 . In FIGS. 3A-3E, the fitted curves are only slightly quadratic and 

appear to be nearly linear in shape. Although, the curves in FIGS. 3A-3E 
are slightly quadratic and nearly linear, the shape of the curves of C v C 2 , 
C 3 , C 4 , and C 5 versus C 0 may generally be linear, quadratic, or any other 
shape as determined empirically or by an analytical method. Thus, while 
the correlations in FIGS. 3A-3E are of quadratic form, any functional 
dependence of C.,, C 2 , C 3 , C 4 , and C 5 on C 0 is within the scope of the 

present invention. 

' 0063 ' Based on the curves depicted in FIGS. 3A-3E, C„ C 2 , C 3 , C 4 , and C 5 are 
each correlated with C 0 and there is thus only one degree of statistical 
freedom among C 0 , C v C 2 , C 3 , C 4 , and C 5 . In using one degree of 



statistical freedom to simulate uncertainties such as process variations by 
Equation (2), the method of the present invention picks C 0 randomly from 
a probability density function (PDF) of C 0 and uses correlations such as 
shown in FIGS. 3A-3E to determine C v C 2 , C 3 , C 4 , and C 5 from the picked 
value of C 0 . A result from employing such a technique is shown in FIG. 4, 
which shows total capacitance versus applied voltage curves, in 
accordance with embodiments of the present invention. In the tests from 
which the curves of FIG. 4 are obtained, the total capacitances of 5 
devices on each of 42 dies were measured from which the data shown in 
FIG. 4 are drawn. This measured data was used to extract area and 
perimeter capacitances. Each die of the 42 dies had five, devices and one 
device was selected from each die for the data shown in FIG. 4. All 5 
devices were located close to one another (i.e., within about 3000 
microns) on the same die. The anode to cathode voltage of the varactor is 
swept from slightly forward bias to relatively large reverse bias but less 
than the onset of avalanche multiplication. 

[0064] 

In FIG. 4, dashed lines 30 represent test data exhibiting characterized by 
process variations, line 32 represents a nominal fit of Equation (2) to the 
test data, and bounding curves 34 and 36 representing the limiting 
envelope of simulated process variations in conjunction with Equation (2) 
by picking C 0 randomly from its PDF and calculating G,, C 2 , C 3 , C 4 , and 
C 5 from the picked value of C 0 through the correlations analogous to those 
shown in FIGS. 3A-3E. FIG. 4 shows that curves 34 and 36 envelope the 
test data in a manner that confirms the accuracy of the method of the 



present invention in using Equation (2) to model the test data. 

[0065] The PDF for C 0 may be determined as follows. Capacitance C 0 at zero 
applied voltage may be measured at a multiplicity of points on a wafer 
and/or on different wafers, resulting in a corresponding multiplicity of 
capacitances C 0 which reflect a variety of process variations. The 

measured multiplicity of capacitances may be fitted to a probability density 
function such, inter alia, as a normal probability distribution by any method 
known to a person of ordinary skill in the art (e.g., by calculating the mean 
of the normal distribution as the arithmetic mean of the measured 
capacitances, and by calculating the standard deviation of the normal 
distribution as the standard deviation of the measured capacitances). 

[0066] Thus a method for developing a model for incorporating process variability 
comprises determining a PDF for C 0 ; and correlating C v C 2 , C 3 , C 4 , and 
C 5 with C 0 as discussed supra in conjunction with FIGS. 3A-3E. 

[0067] A method for using the developed model for randomly selecting a value of 
capacitance C comprises: providing a voltage value V=V", picking C 0 from 
its PDF; calculating C v C 2 , C 3 , C 4 , and C 5 from their correlation with C 0 , 
and substituting C 0 , G,, C 2 , C 3 , C 4 , C 5 and V" into Equation (2) to 
determine C. 

[0068] 

A method for assessing a design that includes capacitance subject to 
process variations comprises: evaluating a performance characteristic of 
the design (e.g., signal transmission speed, signal switching speed, 
operating temperature, power output, power conversion efficiency, the 



logical correctness of a logic circuit, etc.) by simulating the design. 
Simulating the design involves calculating capacitance at nodes (i.e., 
circuit locations) in the design by: providing the voltage V=V(i) at the node 
i, picking a value of C 0 from its PDF; calculating C v C 2 , C 3 , C 4 , and C 5 
from their correlation with C 0 ; and substituting C 0 , C v C 2 , C 3 , C 4 , C 5 and V 
(i) into Equation (2) to determine C. The determined values of C at the 
nodes are utilized to determine the performance characteristic. For 
example, the performance characteristic of switching speed may be a 
function of the values of capacitances C(1), C(2), ... so determined. 

[0069] While the preceding methodology picked C 0 randomly, any of the other 
parameters C v C 2 , C 3 , C 4 , and C 5 could have alternatively been randomly 
picked instead of C 0 . For example, C 1 could have been picked randomly 
from its PDF and the other parameters C 0 , C 2 , C 3 , C 4 , and C 5 would be 
calculated from their correlations with C r 

[0070] 

While the fifth order polynomial of Equation (2) has been utilized for 
implementing the method of the present invention, any equation that 
acceptably models or fits the characteristic (e.g., capacitance) of interest 
may alternatively be utilized. For examples, Equation (1) is suitable for 
calculating capacitance in some applications and may be utilized as 
follows. A probability density function (PDF) may be determined for C 0 , 
and parameters V B and m by measuring correlated with C 0 by test 
measurement analogous to those obtained in conjunction with FIGS. 3A- 
3E, described supra. The capacitance C in Equation (1) may be randomly 
sampled by picking a value of C 0 from its PDF; calculating V B and m from 



their correlation with C 0 ; and substitute the picked value of C 0 and the 
calculated values of V B and m into Equation (1) to determine C. 

[0071] A design that comprises capacitance may be assessed by: evaluating a 
performance characteristic of the design by simulating the design, which 
involves calculating capacitance at nodes (i.e., locations) in the design by: 
providing the voltage V=V(i) at node i, picking a value of C 0 from its PDF; 
calculating V B and m from their correlation with C 0 ; and substituting V B 
and m into Equation (1) to determine C. 

[0072] In another approach, the capacitance of a physical or virtual capacitor may 
be expressed in terms of area capacitance C A and perimeter capacitance 
C p via Equation (3). 

[0073] C = A*C A + P*C p (Equation 3) 

[0074] wherein A is the plate area and P is the plate perimeter of the capacitor, 
wherein A and P are assumed to be known at each location (e.g., in the 
wafer) where C is to be calculated, C A and C p may each be modeled by 

Equation (2) or any other applicable fit equation (e.g, Equation (1) in some 
applications); i.e., 

[0075] C A = C 0A (1 + C 1A V + C 2A V 2 + C 3A V 3 + C^V 4 + C 5A V 5 ) (Equation 4) 
[0076] C p = C op (1 + C 1P V + C 2P V 2 + C 3P V 3 + C^V 4 + C 5p V 5 ) (Equation 5) 
[0077] 

An algorithm of the present invention: determines C A and C p by random 
sampling; and substituting C A and C p into Equation (3) to obtain C. C A and 
C p may be calculated independently or correlatively, depending on 



whether C 0A and C op are sufficiently correlated. A correlation between C 0A 
and C op may be deduced from scatter data of C 0A as a function of C 0P> 
wherein the points of the scatter data may represent pairs of (C 0A , C 0P ) at 

a variety of locations on a wafer and/or several wafers. For a given applied 
voltage V, each pair of (C A , C p ) may be deduced from Equation (3), by 

measuring C for at least two sufficiently close locations on the wafer (e.g., 
on a kerf surrounding a semiconductor chip) at which C A and C p are not 

expected to vary due to process variation. Therefore at zero applied 
voltage V, each pair of (C 0A , C op ) may be deduced from Equation (3) by 

measuring C (at zero volts) for least two sufficiently close locations on the 
wafer at which C 0A and C 0P are not expected to vary due to process 
variation. Thus, (C A> C p ) and (C 0A , C op ) may generally be inferred from 
measurement of total capacitance C for at least two sufficiently close 
locations on the wafer, at a given non-zero applied voltage V and a given 
zero applied voltage V, respectively. 

[0078] 

A correlation coefficient r (e.g., a Pearson correlation coefficient) may be 
calculated from the scatter data of C 0A versus C op . Although the scatter 

data may be represented graphically in the form of a scatter plot, an actual 
scatter plot is not required for calculation the correlation coefficient r, since 
the correlation coefficient r may be calculated directly from the scatter data 
itself. Since the correlation coefficient r is generally in a range of -1 to 1, it 

is useful to define a correlation parameter R = r 2 , which falls in the range 
of 0 to 1 . C 0A and C 0P may be considered as being sufficiently correlated if 

the correlation parameter R is no less than a specified value R MIN (e.g., 



R MIN may be in a range of, inter alia, 0.70 - 0.90). If C 0A and C 0P are 
sufficiently correlated, then C 0A may be regressed upon C op to determine 
an equation that correlates C 0A with C 0P . Thus if C 0A and C 0P are 
sufficiently correlated, then C 0P may be picked randomly from its PDF and 
C 0A may then be calculated from the equation that correlates C 0A with C op . 
Alternatively, C 0A may be picked randomly from its PDF and C 0P may then 
be calculated from the equation that correlates C op with C 0A which is the 
inverse of the equation that correlates C 0A with C 0P . However, if C 0A and 
C op are not sufficiently correlated, then C 0A and C op may be independently 
picked randomly from their respective PDFs. 

[0079] 

The coefficients C 1A , C^, C 3A , C 4A , C 5A of Equation (4) are each 
correlated with C 0A in a manner analogous to the correlation of C v C 2 , C 3 , 
C 4 , and C 5 with C 0 as described supra in conjunction with FIGS. 3A-3E. 
Similarly, the coefficients C 1P , C 2P , C 3P , C 4P , C 5P of Equation (5) are each 
correlated with C 0P in a manner analogous to the correlation of C v C 2 , C 3 , 
C 4 , and C 5 with C 0 as described supra in conjunction with FIGS. 3A-3E. 
After values of C 0A and C op are determined as described supra, the 
remaining sets of coefficients (C 1A , C 2A , C 3A , C 4A , C 5A ) and (C 1P , C 2P , C 3P , 
C 4P , C 5P ) are determined from their correlation with C 0A and C 0P , 
respectively. Then C A and C p may be calculated from Equations (4) and 

(5), respectively, which enables C to be subsequently determined by 
substituting C A and C p into Equation (3). The determined values of C at 

the locations in the design are then utilized to determine the performance 
characteristic. For example, the performance characteristic of switching 



speed may be a function of the values of capacitance C so determined. 

[0080] The preceding discussion of FIGS. 1-4 and Equations 1-5 for the example 
of the capacitance of a varactor will next be extended in the second and 
third segments of the description herein to formulate the scope of the 
present invention in its generality. The second segment of the description 
herein formulates the present invention generally by discussing the 
modeling of a characteristic C that is distributed within a domain, wherein 
C is not expressed in terms of any subcharacteristics of C. The third 
segment of the description herein formulates the present invention 
generally by discussing the modeling of a characteristic C that is 
distributed within a domain, wherein C is expressed in terms of a plurality 
of subcharacteristics of C. 

[0081] Modeling a Characteristic Not Expressed in Terms of Subcharacteristics 

[0082] The present invention discloses a method for statistically modeling a 

characteristic C that is distributed within a domain, wherein C is a function 
of a variable V, and wherein C is not expressed in terms of 
subcharacteristics, as embodied in the example of Equation (2) discussed 
supra. The domain may be a physical domain (e.g., a physical surface or 
volume) such as, inter alia, a semiconductor wafer comprising integrated 
circuits, a geographic neighborhood comprising houses, a galaxy 
comprising stars, etc. For the example of a semiconductor wafer 
comprising integrated circuits, the characteristic C may be capacitance 
and the variable V may be applied voltage, as discussed supra in 
conjunction with Equation (2). The domain may also be a non-physical 



domain such as, inter alia, the domain of stocks on the New York Stock 
Exchange, wherein the characteristic C may be, interaUa, the price-to- 
earnings (PE) ratio of each stock and the variable V may be, inter alia, an 
interest rate benchmark (e.g., the Federal Funds interest rate). 

[0083] The method is next described in conjunction with FIGS. 5-7, in accordance 
with embodiments of the present invention. FIG. 5 is a flow chart for 
providing a model for representing the characteristic C. FIG. 6 is a flow 
chart for randomly sampling C based on the model for C described by 
FIG. 5. FIG. 7 is a flow chart for determining a performance characteristic 
of a design, based on the model of C in FIG. 5 and the method for 
randomly sampling C in FIG. 6. 

[0084] FIG. 5 is a flow chart 40 depicting steps 41-43 which describe providing a 
model for C. Step 41 provides a base equation for modeling C. The base 
equation expresses the characteristic C as a function f of the variable V 
through use of N+1 parameters C 0 , C N . The base equation is of the 

form: 

[0085] C = f(C 0 , C v .... C N , V) (Equation 6) 

[0086] wherein N is at least 1 and parameters C 0 , C 1 C N are subject to 

uncertainty. A specific example of the base equation (6) is Equation (2) 
wherein N=5. Another example of the base equation (6) is Equation (1) 
wherein N=2 (i.e., C,=V B and C 2 =m). 

[0087] g te p 42 provides a probability density function (PDF) for describing the 
probability of occurrence of C 0 in accordance with the uncertainty in C 0 . 



The PDF may be a normal probability distribution or any other form of 
probability distribution. The discussion supra provided an explanation of 
how the PDF of C 0 relating to the capacitance of Equation (2) may be 

determined from an analysis of test data. This methodology for providing 
the PDF of C 0 from test data may be formulated by postulating that the 

domain comprises nodes such that test data is obtained at the nodes. For 
the capacitance example discussed supra, the nodes may each 
correspond to an integrated circuit on a wafer. 

[0088] To determine the PDF of C 0 from testing, test data of C = C 0 (k) may be 
obtained at each node k of K nodes (K at least 2) such that V = V 0 at node 
k, wherein V 0 = 0 volts in the capacitance example discussed supra, 
wherein C 0 (k) is defined as being C 0 at node k, and wherein k has values 

of 1 , 2 K. The PDF of C 0 is derived from the aforementioned test data 

as described supra, i.e., by any method known to a person of ordinary skill 
in the art of fitting a functional form of a PDF to test data. 

[0089] The discussion supra provides an explanation of how the PDF of C 0 may 

be determined from an analysis of test data. However, the scope of the 
present invention includes determining the PDF of C 0 not only from test 

data, but alternatively by any other method such as from an analytical 
calculation or from computer simulation based on an assumed analytical 
model. 

[0090] Step 43 provides subsidiary equations expressing C 1 C N in terms of 

C 0 . An example of subsidiary equations are equations representing the 
dependencies depicted in FIGS. 3A-3E, discussed supra, which have the 



general form of: C x = g x (C x ^) for functions g x such that x has values of 1 , 

2 N. Although in FIGS. 3A-3E specifically, C x is a linear of quadratic 

function g x of C x _ v the parameter C x may generally be any function g x of 
C x1 which accurately represents the functional relationship between C x 
and C x r Note that the equations C x = g x (C x _.,) may be transformed into 
the equivalent form: C x = h x (C 0 ) for x = 1, 2, .... N. 

[0091] To determine the subsidiary equations from testing, test data of C(k) 

versus V may be obtained at each node k of K nodes, wherein K is at least 
2, wherein a "node" is a different location in the region of testable devices 
(e.g., test sites on a wafer), wherein C(k) is defined as C at node k, and 
wherein k has values of 1 , 2, K. The function f in Equation (6) is fitted to 
the test data at each node k to obtain C(k) = f(C 0 (k), C^k), .... C N (k), V), 
wherein C 0 (k), C^k), .... C N (k) respectively denote C 0 , C v C N at node 

k. The subsidiary equations are subsequently derived by utilizing C(k) = f 
(C 0 (k), C,(k) C N (k), V) at each node k of the K nodes. For example, a 

curve of a given analytical form (e.g., linear, quadratic, exponential, etc.) 
for representing C 1 versus C 0 can be fitted to the K points (C Q ^),C^)), 
(C 0 (2),C 1 (2)), .... (CqW.C^K)) 

[0092] FIG. 6 is a flow chart 45 depicting steps 40 and 46-49 for randomly 
sampling C. Step 40 is the flow chart 40 of FIG. 5 which provides the 
model for C, described supra. Given the model for C provided by step 40, 
steps 46-49 are executed when a value of C is to be randomly sampled. 

[0093] Step 46 provides a value V" of V as input to the process of randomly 
sampling C. 



[0094] Step 47 picks a random value C 0R of C 0 from the PDF of C 0 , wherein the 

PDF was provided in step 42 of FIG. 5 as described supra. Given a PDF 
of C 0 , any method for picking C 0R from the PDF may be utilized (e.g., 
picking C 0R by utilizing the cumulative probability distribution associated 
with the PDF as is known in the art). 

[0095] Step 48 computes values C 1R , .... C NR of C.,, C N , respectively, by 

substituting C 0R into the subsidiary equations that were provided in step 
43 of FIG. 5, described supra. 

[0096] Step 49 calculates C by substituting C 0R , C 1R , .... C NR and V" into the base 
equation provided in step 41 of FIG. 5. 

[0097] FIG. 7 is a flow chart 50 depicting steps 40 and 51-52 for determining a 
performance characteristic of a design, based on the model of C in FIG. 5 
and the method for randomly sampling C in FIG. 6. The design comprises 
I nodes in the domain in which C is distributed, wherein I is at least 2. 
Each node i of the I nodes has a value C(i) of the characteristic C, wherein 
i has values of 1 , 2, .... I.. Step 40 is the flow chart 40 of FIG. 5 which 
provides the model for C, described supra. 

[0098] Given the model for C provided by step 40, step 51 randomly selects a 
value of C(i) of C at each node i of the I nodes. At each node i, step 51 is 
implemented by executing steps 46-49 of FIG. 6, wherein V" in step 46 is 
V(i). V(i) is defined as the value of V at node i, which is input to the 
calculation of C(i) as discussed supra in conjunction with step 46 of FIG. 
6. 



[0099] Given the values of C(1), C(2) C(l) randomly selected in step 51 , step 

52 determines the performance characteristic, by utilizing the randomly 
selected values of C(1 ), C(2) C(l). 

[0100] Modeling a Characteristic Expressed in Terms of Subcharacteristics 

[0101] The present invention discloses a method for statistically modeling a 

characteristic C that is distributed within a domain, wherein C is a function 
of a variable V, and wherein C is expressed in terms of J 
subcharacteristics S 1f S 2 S Jt as embodied in the example of Equation 

(3) discussed supra. The domain may be a physical domain(e.g., a 
physical surface or volume) such as, inter alia, a semiconductor wafer 
comprising integrated circuits, a geographic neighborhood comprising 
houses, a galaxy comprising stars, etc. For the example of a 
semiconductor wafer comprising integrated circuits, the characteristic C 
may be capacitance, the subcharacteristics of C may be area capacitance 
C A and perimeter capacitance C p (i.e., J=2), and the variable V may be 

applied voltage, as discussed supra in conjunction with Equation (3). The 
domain may also be a non-physical domain such as, inter alia, the domain 
of stocks on the New York Stock Exchange, wherein the characteristic C 
may be, interaWa, the price-to-earnings (PE) ratio of each stock, the 
subcharacteristics may be stocks paying dividends and stocks not paying 
dividends (i.e., J=2), and the variable V may be, inter alia, an interest rate 
benchmark (e.g., the Federal Funds interest rate). 

[0102] The method is next described in conjunction with FIGS. 8-10, in 

accordance with embodiments of the present invention. FIG. 8 is a flow 



chart for providing a model for representing the characteristic C in terms of 
its J subcharacteristics S v S 2 Sj. FIG. 9 is a flow chart for randomly 

sampling C based on the model for C and its subcharacteristics described 
by FIG. 8. FIG. 10 is a flow chart for determining a performance 
characteristic of a design, based on the model of C and its 
subcharacteristics in FIG. 8 and the method for randomly sampling C in 
FIG. 9. 

[0103] FIG. 8 is a flow chart 60 depicting steps 61-65 which describe providing a 
model for C. Step 61 provides a combination equation for modeling C as a 
function F its J subcharacteristics S r S 2 Sj. The combination equation 

is of the general form: 
[0104] C = F(S 1( S 2> .... Sj) (Equation 7) 

[0105] as illustrated by the example of Equation (3) in which C is a linear 
combination of C A and C p (i.e., J=2, S,,=C A , S 2 =C p ). 

[0106] Step 62 provides base equations which express Sj as a function f. of a 
variable V through use of N+1 parameters S j0 , S j1f S jN . The base 
equations are of the form: 

[0107] ^ = fj(Sj 0> S jv .... S jN , V) (Equation 8) 

[0108] wherein N is at least 1. The parameters S j0 , S^, .... S jN are subject to 

uncertainty, and the subscript j has values of 1 , 2, J. For a fixed value 
of j, the parameters Sj 0 , S^, .... S jN in base equations (8) are analogous to 
parameters C 0 , C 1 C N in Equation (6). 



[0109] In Equation (8), l } may have various functional characteristics. A first 
example of a functional characteristic for fj is that fj may have a same 
functional form of V for each subcharacteristic Sj of the J 
subcharacteristics (e.g., a fifth order polynomial in V for each S.). A 
second example of a functional characteristic for fj is that f. may be 
constant with respect to any variation in V for each subcharacteristic S. of 
the J subcharacteristics (i.e., fj does not depend on V). A third example of 
a functional characteristic for fj is that fj may vary with respect to a 
variation in V for each subcharacteristic Sj of the J subcharacteristics (i.e., 
each fj varies with V). 

[01 10] Step 63 provides at least one probability density function (PDF) selected 

from PDF r PDF 2 PDFj. PDF n describes the probability of occurrence 

of S n0 for n=1, 2, J. 

[011 1] PDF r PDF 2 , .... PDFj may each be a normal probability distribution or any 
other form of probability distribution. At least one PDF n is provided and 
PDF 1 is always provided. Thus for n>0, none, some, or all of the PDF n 
may be provided. If provided, PDF n may be determined by any method 
described supra for determining the PDF of C 0 , including test methods, 
analytical calculation, or computer simulation based on an assumed 
analytical model. 

[0112] As stated supra, PDF 1 is always provided. PDF n (for n>0) is not provided if 
S n0 is sufficiently correlated with S 10 . S n0 is sufficiently correlated with S 10 
if a correlation parameter R n between S n0 and S 10 is no less than a 
specified minimum correlation parameter R MIN . The correlation parameter 



R n is the square of a correlation coefficient r n between S n0 and S 10 The 
correlation coefficient r n may be determined from scatter data of S n0 
versus S 10 as described supra. 

[01 13] For each PDF n not provided, step 64 provides an auxiliary equation E n 
expressing S n0 in terms of S 10 which is analogous to the aforementioned 
relating of C 0A to C 0P to determine an equation correlating C 0A with C 0P , as 
discussed supra. The auxiliary equation E n may be derived from the 
scatter data of S n0 versus S 10 . 

[01 14] Step 65 provides subsidiary equations expressing S )V .... S jN in terms of 
S j0 for each subcharacteristic Sj of the J subcharacteristics, which for a 
fixed j are analogous to the subsidiary equations expressing C v .... C N in 
terms of C 0 as described supra in conjunction with FIGS. 3A-3E. Thus the 
subsidiary equations for calculating S^, .... S jN in terms of S j0 (for fixed j) 

may be derived by any of the methods described supra tor deriving the 
subsidiary equations for calculating C v C N in terms of C 0 . 

[01 15] FIG. 9 is a flow chart 70 depicting steps 60 and 71-75 for randomly 
sampling C. Step 60 is the flow chart 60 of FIG. 8 which provides the 
model for C in terms of its J subcharacteristics S v S 2 Sj, described 

supra. Given the model for C provided by step 60, steps 71-75 are 
executed when a value of C is to be randomly sampled. 

[01 16] Step 71 provides a value V" of V as input to the process of randomly 
sampling C. 

[01 17] For each poF n provided in step 63 of FIG. 8 (including PDF.,), step 72 



picks a random value of S n0R from PDF n . For each PDF n not so provided, 
step 72 calculates S n0R by substituting S 10R into the auxiliary equation E n 
provided in step 64 of FIG. 8. 

[01 18] Step 73 computes values S j1R , .... S jNR of S jN , respectively, by 

substituting S j0R into the subsidiary equations provided by step 65 of FIG. 
8 for each subcharacteristic Sj of the J subcharacteristics. 

[01 19] Step 74 calculates Sj by substituting S j0R , S j1R , .... S jNR and V" into the 

base equations provided by step 62 of FIG. 8 for each subcharacteristic Sj 
of the J subcharacteristics 

[0120] Step 75 calculates C by substituting S v S 2 , .... Sj into the combination 
equation provided by step 61 of FIG. 8. 

[0121] FIG. 10 is a flow chart 77 depicting steps 60 and 78-79 for determining a 
performance characteristic of a design, based on the model of C in FIG: 8 
and the method for randomly sampling C in FIG. 9. The design comprises 
I nodes in the domain in which C is distributed, wherein I is at least 2. 
Each node i of the I nodes has a value C(i) of the characteristic C, wherein 

i has values of 1 , 2 I.. Step 60 is the flow chart 60 of FIG. 8 which 

provides the model for C in terms of its J subcharacteristics S v S 2 Sj, 

described supra. 

[0122] Gjven the model for C in terms of its J subcharacteristics S r S 2 , Sj 

provided by step 60, step 78 randomly selects a value of C(i) of C at each 
node i of the I nodes. At each node i, step 78 is implemented by executing 
steps 71-75 of FIG. 9, wherein V" in step 71 is V(i). V(i) is defined as the 



value of V at node i, which is input to the calculation of C(i) as discussed 
supra in conjunction with step 71 of FIG. 9. 

[0123] Given the values of C(1), C(2),..., C(l) randomly selected in step 78, step 
79 determines the performance characteristic, by utilizing the randomly 
selected values of C(1), C(2), .... C(l). 

[0124] While the description of embodiments of the present invention has 
focused on modeling a characteristic C or subcharacteristic Sj that 

depends on a single independent variable V, the scope of the present 
invention also includes modeling C and Sj each as a function of a plurality 

of independent variables V 1f V 2 wherein the base Equations (6) and 

(8) for C and Sj would be respectively generalized to Equations (9) and 
(10): 

[0125] C = f(C 0 , C 1 C N ,V V V 2 , ...)(Equation 9) 

[0126] S, = fj(Sj 0 , S jv .... S jN , V v V 2 , ..) (Equation 10) 

[0127] -r ne dependence of C and/or Sj on V.,, V 2 , ... may be individually tailored 
for each of V r V 2 , .. and may be cast in any applicable functional form for 

each independent variable such as, inter alia, the polynomial functional 
form of Equation (2) and/or Equations (4)-(5). The parameters C 0 , C v ... 
C N would be coupled to the variables \i v V 2 , ... in a manner that reflects 

the functional forms chosen. For example for the case of fitting C to a 
polynomial in each independent variable of two independent variables V 1 
and V 2 , some of the parameters of C 0 , C v ... may appear as coefficients of 
powers of V 1 whereas other parameters of C 0 , C.,, ... may appear as 



coefficients of powers of V 2 . Similar to the single independent variable V 
embodiment when picking a random value of C, the parameter C 0 for the 
multiple independent variable V v V 2 , ... embodiment would be picked from 
its PDF and parameters C 1( .... C N would be determined based on their 
relationship to C 0 as expressed in subsidiary equations in accordance with 
the description supra thereof. 

[0128] Computer System 

[0129] 

FIG. 1 1 illustrates a computer system 90 used for implementing algorithms 
that model and sample a characteristic C (and assess a design comprising 
C), in accordance with embodiments of the present invention. The 
computer system 90 comprises a processor 91 , an input device 92 
coupled to the processor 91 , an output device 93 coupled to the processor 
91 , and memory devices 94 and 95 each coupled to the processor 91 . The 
input device 92 may be, inter alia, a keyboard, a mouse, etc. The output 
device 93 may be, inter alia, a printer, a plotter, a computer screen, a 
magnetic tape, a removable hard disk, a floppy disk, etc. The memory 
devices 94 and 95 may be, inter alia, a hard disk, a floppy disk, a 
magnetic tape, an optical storage such as a compact disc (CD) or a digital 
video disc (DVD), a dynamic random access memory (DRAM), a read- 
only memory (ROM), etc. The memory device 95 includes a computer 
code 97. The computer code 97 includes one or more algorithms that 
model and sample the characteristic C (and assess a design comprising 
C), in accordance with embodiments of the present invention. The 
processor 91 executes the computer code 97. The memory device 94 



includes input data 96. The input data 96 includes input required by the 
computer code 97. The output device 93 displays output from the 
computer code 97. Either or both memory devices 94 and 95 (or one or 
more additional memory devices not shown in FIG. 11) may be used as a 
computer usable medium (or a computer readable medium or a program 
storage device) having a computer readable program code embodied 
therein and/or having other data stored therein, wherein the computer 
readable program code comprises the computer code 97. Generally, a 
computer program product (or, alternatively, an article of manufacture) of 
the computer system 90 may comprise said computer usable medium (or 
said program storage device). 

[0130] While FIG. 1 1 shows the computer system 90 as a particular configuration 
of hardware and software, any configuration of hardware and software, as 
would be known to a person of ordinary skill in the art, may be utilized for 
the purposes stated supra in conjunction with the particular computer 
system 90 of FIG. 1 1. For example, the memory devices 94 and 95 may 
be portions of a single memory device rather than separate memory 
devices. 

[0131] While embodiments of the present invention have been described herein 
for purposes of illustration, many modifications and changes will become 
apparent to those skilled in the art. Accordingly, the appended claims are 
intended to encompass all such modifications and changes as fall within 
the true spirit and scope of this invention. 



